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Abstract: Modifications of the initial-state of the inflaton field can induce a departure 
from Gaussianity and leave a testable imprint on the higher order correlations of the CMB 
and large scale structures in the Universe. We focus on the bispectrum statistics of the 
primordial curvature perturbation and its projection on the CMB. For a canonical single- 
field action the three-point correlator enhancement is localized, maximizing in the collinear 
limit, corresponding to enfolded or squashed triangles in comoving momentum space. We 
show that the available local and equilateral template are very insensitive to this localized 
enhancement and do not generate noteworthy constraints on initial-state modifications. 
On the other hand, when considering the addition of a dimension 8 higher order derivative 
term, we find a dominant rapidly oscillating contribution, which had previously been over- 
looked and whose significantly enhanced amplitude is independent of the triangle under 
consideration. Nevertheless, the oscillatory nature of (the sign of) the correlation func- 
tion implies the signal is nearly orthogonal to currently available observational templates, 
strongly reducing the sensitivity to the enhancement. Constraints on departures from the 
standard Bunch-Davies vacuum state can be derived, but also depend on the next-to- 
leading terms. We emphasize that the construction and application of especially adapted 
templates could lead to CMB bispectrum constraints on modified initial states already 
competing with those derived from the power spectrum. 
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1. Introduction 

It has been known for a while now [1,2] that the potential presence of non-Gaussian signa- 
tures in the CMB is a powerful probe of the physics of inflation and beyond. Computations 
of the primordial bispectra [3,4, 7-12, 14] (and later trispectra [2, 18]) have shown that dif- 
ferent models of inflation can produce rather unique features, which would allow, when 
detected, to discriminate between them. For the bispectrum, the distinction between mod- 
els relies on two features 1) the overall amplitude of the non-Gaussian signal and 2) the 
detailed dependence on the comoving momenta. Obviously, when the overall amplitude of 
the signal is low, the second feature will be much harder to observe. Observational lim- 
itations due to foreground contamination [20,21] and cosmic variance limit the detection 
of non-Gaussianity in the CMB temperature and polarization spectrum [22, 24] . For that 
reason one can already conclude that non-Gaussianity should be observably absent if a 
single, slowly rolling, scalar field is responsible for inflation [25,26]. 

Even if a model predicts a detectable non-Gaussian amplitude, it will remain a chal- 
lenge to measure the actual momentum dependence, since the inferred constraints on the 
level of non-Gaussianity [33-37] are based on a sum over all modes of a pre-assumed momen- 
tum dependence 1 . Such dependencies are known as 'local' (or 'squeezed') and 'equilateral' 

1 Other methods such as Minkowski Functionals (see [31] for theory and [37] and references therein for 
observational results) and a Wavelet approach ( [32] and subsequent papers) exist which typically do not 
rely on a pre-assumed momentum dependence. Here however we refer to the approach initiated in [22] and 
further developed in [39-46], which seems to give the most consistent and stringent constraints [37] so far. 
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template, which correspond to particular shapes that maximize in some 'extreme' triangle 
configuration in momentum space. The possibility to distinguish between different theo- 
retical models producing a sizable non-Gaussian amplitude relies on the fact that in the 
models considered so far the produced non-Gaussianities are well approximated by one of 
these templates. For example, it has been shown [14,17] that non-canonical kinetic terms 
and higher derivative contributions to the inflaton potential can produce significant levels 
of non-Gaussianity of the equilateral type if the speed of sound in these models is much 
smaller than the speed of light, which can be realized in certain brane inflation scenar- 
ios [13,15]. Local shape non-Gaussianities were the first type to be considered [5,22,38] 
and are a direct consequence of the nonlinear relation between the inflaton fluctuations 
and the curvature perturbations that couple to matter and radiation. In [25, 26] it was 
shown that the amplitude of local type non-Gaussianities in single-field slow-roll inflation 
is proportional to the slow-roll parameter e 2 , which is very small by construction. The 
amplitude of local type non-Gaussianities should therefore be undetectably small if single- 
field slow-roll is responsible for inflation. In contrast, large local non-Gaussianities can be 
generated in curvaton models [7], where the curvature perturbation ( can evolve outside 
the horizon, or inflationary models with multiple scalar fields. Models of new ekpyrosis [11] 
and the recently proposed contracting models with an increasing speed of sound [17] , in 
which a bouncing universe is replacing inflation, also yield large non-Gaussianities of the 
local type. 

In this paper we will focus on non-Gaussian features arising from an arbitrary initial- 
state modification. This type of non-Gaussianity has been discussed in [49-54], here we 
provide a more detailed analysis on their detectability. In the language of boundary effective 
field theory [29] one can generally divide the contributions into two parts; non-Gaussianities 
coming directly from the initial-state boundary (which are absent when considering Gaus- 
sian initial-state modifications), and 'bulk' non-Gaussianities generated by the presence of 
(interacting) particles in the modified initial-state [55] . In the boundary effective field the- 
ory formalism it has been shown that the leading non-Gaussian initial-state modification 
is of the local type [51,52]. However, the 'bulk' non-Gaussianities generated by the non- 
zero Bogoliubov coefficient seem to have a unique momentum dependence, which is very 
different from that of the local and equilateral types [14, 54, 55] . For example, for a canon- 
ical single-field inflaton action, in momentum space the non-Gaussian signal produced by 
a, possibly Gaussian, modification of the initial-state maximizes for triangles where two 
momentum vectors are collinear, i.e. when the magnitude of one of the comoving momenta 
equals the sum of the other two: k p = k q + k r with p / q / r (a squashed, flattened or 
enfolded triangle), and are known as collinear or enfolded type non-Gaussianities. 

Thus far only the local and equilateral type of non-Gaussianities have been constrained 
by the data [35,37], although recently a strong case has been made for a more general 
set-up [57]. There are essentially two reasons for this. First of all, the realization that 
initial-state modifications give rise to a unique non-Gaussian shape, that might even be 
detectable due to subtle enhancements, is rather recent and its theoretical motivation might 

2 This is the semi-classical, tree level, result. The effects of quantum loop corrections have been studied 
in [6]. 
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be considered less compelling. Putting aside plausible theoretical concerns associated with 
modifications of the vacuum state and instead taking a phenomenological point of view, 
deviations from the standard Bunch-Davies state are tightly constrained [23] because of 
their unique oscillatory signatures in the 2-point power spectrum (see [30] and references 
therein). As emphasized in [54] the bispectrum (and possibly even higher n-point func- 
tions) might be as good, or even better, in constraining initial-state modifications. Clearly, 
with the expected future improvements in detecting primordial non-Gaussian signals, it 
is worthwhile to look for the presence of enfolded type or, as we will argue in this paper, 
oscillatory non-Gaussianities in the CMB data to constrain initial-state modifications. The 
second more pragmatic reason why enfolded type non-Gaussianities have not been com- 
pared to the data yet is that in analyzing the data computational limitations demand that 
the momentum dependence is factorizable. Generic 3-point correlators are not factorizable, 
so one resorts to constructing a factorizable template that approximates the actual theo- 
retical bispectrum, maximizing in the appropriate 'extreme' triangle. It is this factorizable 
template that is then compared to the data. Such templates have been constructed for 
the local [2] and equilateral shapes [34], but has not yet been constructed for the type 
of non-Gaussianities predicted by initial-state modifications, which are typically expected 
to extremize in an enfolded (collinear, squashed or flattened) triangle. The goal of this 
paper is two-fold: to present a detailed analysis on the detectability of non-Gaussianities 
produced by initial-state modifications using currently available templates, both with and 
without higher derivative corrections, and secondly to determine how much improvement 
can theoretically be gained by using more optimal templates. 

The paper is organized as as follows. In Section 2 we review the standard analytical 
tools to study non-Gaussianity, in particular the computation of the 3-point correlation 
function in momentum space and its relation to different triangular shapes. In Section 3 
we will present a detailed analysis of the 3-dimensional bispectrum from initial-state mod- 
ifications in the single- field slow-roll inflationary scenario. In Section 4 we analyze the 
case of modified initial-state non-Gaussianities in the presence of a dimension 8 higher 
order derivative term in the Lagrangian. In Section 6 we discuss the results of the CMB 
bispectrum computation and finally we present our conclusions in Section 7. 

2. Three-dimensional bispectrum preliminaries 

In this section we will briefly review the standard tools for analyzing non-Gaussianities 
as first described in [40]. In the next sections we will apply these tools to the case of 
initial-state modifications. Let us start considering the primordial spectrum of curvature 
perturbations generated by the inflaton. In three-dimensional comoving momentum space 
a generic three-point correlator of the curvature perturbation is a function of the three 
comoving momenta fci, &2 and k$, which in 3 dimensions corresponds to a total of 9 
parameters. Translational invariance forces the three-point function to conserve momentum 
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which fixes one of the momenta, reducing the number of free parameters from 9 to 6. Ro- 
tational invariance allows one to pick a 2-dimensional plane defined by the remaining two 
momenta and adjust the axes such that one of the momenta is along one of the axes of the 
plane. This fixes another 2 + 1 = 3 parameters, leaving only 3 variables to parametrize 
the three-point correlator. These can be identified with two angles and the overall scale of 
the triangle formed by k%, &2> k 3 . Since the primordial power spectrum is approximately 
scale invariant, we expect the correlator to be a homogeneous function F of degree —6 in 
comoving momentum space, i.e. F(\k±, Xk 3 ) = X~ 6 F(ki, k,2, k 3 ). So (approximate) 
scale invariance fixes the dependency of the three-point correlator on the scale of the trian- 
gle, further reducing the number of free parameters to the 2 angles. Instead of writing the 
function F in terms of these angles, it is most convenient to consider the two independent 
ratios given by the magnitudes of the comoving momenta x<i = /C2/&1 and x 3 = k 3 /k±. In 
order to determine the relevant X2, x 3 domain, one assumes k\ > k<i > ki, giving X2 < 1 
and x 3 < 1, and then uses the triangle constraint to find that 1 — X2 < x 3 < 1, identifying 
the top-right triangle in X2, x 3 space (see e.g. figure 4). Since the distributions are sym- 
metric in X2 an x 3 , one could further reduce the domain by half only considering x 3 > X2- 
Hitherto, unlike the power spectrum, which only depends on the 'reciprocal distance' be- 
tween two-points, the bispectrum F depends on two variables, typically represented by the 
ratios of the magnitudes of the comoving momenta F = F(x2, £3). 

To measure the overall amplitude A in Eq. (2.1), one assumes a particular theoretical 
template shape function F(x2,x 3 ), sums over all triangles and then normalizes appropri- 
ately, taking into account the variance of a given mode in Fourier space. This procedure 
leads to the following estimator of the non-Gaussian amplitude A 



£^ 2 (fci,fc2,fc 3 )/(<<< 

Here the represent the variances of the different modes and the sum runs over all triangles 
in momentum space. The above estimator naturally defines a scalar product between two 
distributions Fx and Fy as [40] 

rp jp ST^ F x(k 1 ,k 2 ,k 3 )F Y (k 1 , k 2 ,k 3 ) 

F X .F Y = ^ -g-5-g . (2.3) 

r fcl /C2 /C3 

This scalar product allows us to quantitatively verify how well a particular template distri- 
bution, say Fx, can be used to constrain a theoretical signal described by the distribution 
Fy ■ In terms of the (reduced set of) parameters X2, X3 the sum over triangles can be written 
as an integral with an appropriate measure equal to x^x^ 

Fx-Fyoc J dx 2 dx 3 Fx(x2,x 3 )Fy(x2,x 3 )x2xl. (2.4) 

To derive optimal constraints using a template Fx one would like the scalar product, or 
the overlap, to be as large as possible. Using the scalar product one can construct a 
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normalization independent 'cosine' between two distributions 

<2 - 5) 

which is close to 1 for shapes that are very similar and considerably smaller than 1 for 
shapes that are very distinct. It follows that optimal constraints can be obtained only 
if the (factorizable) templates, which are used to analyze the data, have a cosine close 
to 1 with the theoretically predicted non-Gaussian signal. Nevertheless, for non-optimal 
templates Fx one can still derive constraints on a theoretically predicted non-Gaussian 
signal Fy provided one introduces the so-called fudge factor Ap, defined as [40] 

The fudge factor allows to deduce the relevant constraints for different theoretical pre- 
dictions Fy using the results inferred from the data analysis of a particular template 
distribution Fx- In such a case the constraint on the amplitude of the type Fx will be 
degraded by a factor 1/Ap, thus the smaller the scalar product between Fy and Fx the 
weaker the constraints on the Fy type non-Gaussianities using the Fx template. Looking 
at Eq. (2.6) it should be clear that optimal constraints can be achieved by maximizing the 
cosine between the template and the theoretical prediction. The other contribution to the 
fudge factor has to do with some conventional choice of normalization for the template 
and the theoretical distribution involved and can be adapted accordingly. We will apply 
these techniques to obtain constraints on non-Gaussianities predicted by modified initial- 
states using the latest results on local and equilateral type non-Gaussianities, and to derive 
what can (theoretically) be gained by analyzing the data with an more optimal (enfolded) 
template. 

Let us briefly discuss the normalization conventions for the non-Gaussian amplitudes, 
which are important for a correct interpretation and comparison of the results obtained 
for different distributions. To compare the local and equilateral template one typically 
equates the distributions in the equilateral triangle ki = ki = kj, [40]. We will follow this 
convention, which allows us to directly use the constraints from the CMB for the local and 
equilateral non-Gaussian amplitudes. To be explicit, for the local template distribution 
the standard definition of the /^l^ parameter, starting from the general three-point func- 
tion in Eq. (2.1), is related to the amplitude A of the three-point function of curvature 
perturbations in the following way 3 

A = (2^) 4 (-^;°L Cal ) (2-7) 

where A$ = g^^^ is the amplitude of the two-point power spectrum, which has been 
observed to be approximately equal to 10~ 10 , M p is the reduced Planck mass and e = 

3 Our sign convention for the non-Gaussian amplitude follows [26,40], which is different from that used 
in [39]. 
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\Mp (yY^j is the first slow-roll parameter. In the above expression for the amplitude A 
we included the overall k% scaling dependence, implying that the local shape i^ local can be 
identified as the following function of the reduced number of variables x 2 , x 3 

F^(x 2 ,x 3 ) = 2(^3 + ^3 +^3) , (2-8) 

\ x 2 x 3 X 2 X 3/ 

For non-Gaussianities of the equilateral type it was shown in [34, 40] that these are well 
approximated by the following shape function 



F e ^(x 2 ,x 3 ) = 6 



11 1 2/1 

3 3 3 3 2 2 \ 2 3 perm 

2 *^3 *^'2"^'3 •^'2*^' 3 \ 2*^*3 



(2.9) 



where the normalization has been fixed such that the local and equilateral template shape 
functions both equal 6 in the equilateral limit x 2 = x 3 = 1. Comparing to the local template 
definition of the non-Gaussian amplitude /nlj this then suggests a similar definition of /nl" 1 

A = (2-) 4 ("'/Mr 1 ) (2-10) 

A crucial property of the local and equilateral template is that they are factorized in their 
comoving momentum dependence. This allows for a drastic (and necessary) reduction in 
the computational time needed to compare the template distributions to the CMB data, 
yielding constraints on the parameters f^i^ an d /nl • The analysis of the WMAP-5 year 
data for local and equilateral non-Gaussianities gives the following limits, [37] 

-9 < /i? L cal < 111 
-151 < /£ q L uil < 253, (2.11) 

which we will use in Section 6. It is worth stressing that the non-Gaussian amplitude 
/nl is not uniquely defined, it depends on a specific choice for the shape function Fx, 
which is equivalent to fixing the integrated norm \Fx\ = V F x ■ Fx ■ It is the combination 
/nl I Fx I that is independent of a particular normalization scheme and which measures the 
(integrated) non-Gaussian amplitude. Obviously any choice will do, as long as one properly 
takes into account the corresponding norm | Fx \ when for example deducing constraints on 
the non-Gaussian amplitude f^ L from the equilateral and local template results. 

3. Modified initial-state non-Gaussianities 

Theoretically predicted three-point functions, evaluated in the regular Bunch-Davies vac- 
uum state, describe a non-Gaussian signal of either local or equilateral type, depending on 
whether higher derivative corrections play a significant role in the inflationary evolution. 
If this is the case, as in DBI models of inflation [15, 16], then the dominant contribution 
is of the equilateral type and can be large enough to be detectable in the near future. 
The existence of different shapes can be nicely understood in terms of the nonlinear origin 
of the non-Gaussian signal. For the local shape it is the nonlinear relation between the 
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inflaton and the curvature perturbation on super-horizon scales that produces the maximal 
effect, whereas in the DBI case nonlinear effects in the inffaton sector are most relevant 
and maximize when all momenta cross the horizon. 

As was shown in [53, 54] non-Gaussian effects can also be generated by dropping the 
assumption that the vacuum state is Bunch-Davies. To fundamentally address the vacuum 
state ambiguity one would first need a full understanding of physics at the highest energy 
scales, where the description in terms of a free inflaton field is expected to break down, as 
well as the physics preceding inflation. Nevertheless departures from the free Bunch-Davies 
state can be studied on a phenomenological basis and it seems worthwhile to use observa- 
tions to constrain the possibilities. The two-point power spectrum already provides strong 
constraints on the initial-state, which has to be close to Bunch-Davies [23]. Interestingly 
though, according to [54], three- and higher n-point functions might be very constraining as 
well, mainly due to subtle enhancement effects, which increase the non-Gaussian amplitude 
in collinear or enfolded triangles. In this section we will focus on the simplest case, with 
a three-point correlation function derived in the general context of slow-roll inflation, but 
evaluated in a vacuum state different from standard Bunch-Davies, as parametrized by an 
undetermined Bogoliubov parameter (3k- The leading non-Gaussian contribution due to 
the appearance of a negative frequency mode is essentially obtained by swapping the sign 
of one of the comoving momenta in the slow-roll inflation result [26]. In appendix B we 
confirm the result first derived in [54] that the correction to the three-point correlation 
function is given by 

3 3 

(3-1) 

In the above expression kj = Yli ki — 2kj and rjo represents the initial conformal time which 
has to be introduced to ensure that the non-Gaussian effects can consistently be calculated 
using an effective field theory description valid below some physical cut-off momentum 
scale M [54], i.e. non-Gaussianities generated by initial-state modifications are sensitive 
to the details of the ultraviolet complete theory and the physical cut-off scale M is intro- 
duced to parameterize our ignorance. This suggests that the initial time rjo should be a 
function of the comoving momentum k, allowing the combination \kr]o(k)\ = M/H 3> 1 
to be a large fixed number independent of k. This prescription treats all comoving mo- 
menta equivalently, tracing back different comoving momenta from the time their physical 
momentum equals the cut-off scale M, preserving scale invariance. Instead considering rjo 
to be some fixed initial time would immediately result in a breaking of scale invariance 
because different comoving momenta would receive contributions from a different range of 
physical momentum scales in the time integral involving the interaction Hamiltonian (see 
appendix B). 

This is very reminiscent of the distinction between the New Physics Hypersurface 
(NPH) [27, 28] and Boundary Effective Field Theory (BEFT) [29, 30] proposals to model 
initial-state modifications. In the latter case one fixes an initial time where one calculates 
corrections to the usual Bunch-Davies initial condition using boundary effective field theory. 
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The result is a Bogoliubov parameter (3k depending on k, resulting in an explicit breaking 
of scale-invariance in the two-point power spectrum. In the NPH scenario one traces every 
momentum mode back to some large physical cut-off scale M and imposes the standard 
flat space vacuum state (corresponding to positive frequency modes only), mode by mode, 
resulting in a prediction for /3& that is independent of k, which only gives rise to a small 
departure from scale-invariance after taking into account the slow-roll evolution of the 
Hubble parameter. Note that the NPH vacuum state proposal is not grounded (yet) in 
some effective field theory scheme that can systematically be applied to calculate quantum 
state corrections, as opposed to the BEFT approach. Both proposals have their problems 
and, as emphasized before, should at this stage and for our purposes here be considered as 
purely phenomenological models distinguished primarily by their consequences for scale- 
invariance. 

For the bispectrum, considering an initial time 770 independent of k (BEFT) or a k- 
independent combination \krjo(k)\ = M/H (NPH) immediately results in a breaking of 
scale-invariance, even independent of any specific ^-dependent prediction for (3k- Because 
the tools to analyze non-Gaussian shapes introduced in the previous section crucially rely 
on scale-invariance of the bispectrum, we will assume that the initial time rjo(k) depends 
on the comoving momentum such that \kr]o\ = M/H is /c-independent, in the spirit of the 
scale-invariant NPH scenario. It would be interesting to relax this assumption and apply 
more general techniques, for instance those recently developed in [57], to analyze the scale- 
dependent non-Gaussianities that arise in a BEFT approach of initial-state modifications. 

Looking at the dependence on comoving momenta, we can see that the three-point 
correlator Eq. (3.1) maximizes in the collinear or enfolded triangle defined by k\/k\ = 
vr / 1 A; 1 770 1 ~ 0, and with an amplitude proportional to k\r]Q. Here k\ is assumed to be the 
largest comoving momentum in the triangle, whose overall power-law dependence mani- 
fests the expected scale invariance of the three-point function. Similarly to the local and 
equilateral shapes, this enfolded type of non-Gaussianity can be associated to a dominant 
source of nonlinearities, in this specific case this is related to the unavoidable presence 
of interacting particles in the modified initial-state at sub-horizon scales. Based on the 
|£;i?7o| = M/H 3> 1 enhancement of the non-Gaussian signal a rough order of magnitude 
estimate of the observational constraints on modified initial-state non-Gaussianities was 
given in [54]. However, their estimate was inferred by considering the maximum signal in 
the enfolded limit, and directly compared to existing bounds on the local non-Gaussian 
amplitude. In contrast, a full analysis of the sensitivity of current non-Gaussian constraints 
on departures from the Bunch-Davies vacuum must involve integrating over all triangles 
and crucially relies on the scalar product between the theoretical template prediction and 
the different observational template distributions. We will address this issue throughout 
the rest of this paper. 

In order to proceed and calculate the scalar product, cosine and fudge factor, we need 
to determine the dominant contribution to the shape function and identify the correspond- 
ing non-Gaussian amplitude. Starting from Eq. (3.1) we identify the relevant comoving 
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momentum dependent part as 

pmodin/, k h\- 1 f 1 ~ COS Mfa + fc 2 - fc 3 )] 1 - COs[7/ (fci - k 2 + fc 3 )] 

1 2 ' 3J fciMsl Jfe§(*i + ^(fci-fca + fcs) 



1 - cos[77o(— h + k 2 + k : 



3 J 



fcf(-fci + A; 2 + fc 3 ) 



(3.2) 



which by having scaled out the standard 6 dependence leads to the corresponding defi- 
nition for the amplitude A, 

A = (2^(3e\(3\)^, 



where we have replaced TZe(Pk) with the absolute value By comparing Eq. (3.3) to the 

f cnf 
'NL 



amplitudes of the local and equilateral templates, a standard definition of /S? suggests 



that — 5e|/3|. Consequently, without any enhancement from a large fudge factor this 
non-Gaussian amplitude is obviously undetectable, since it is suppressed by both the slow- 
roll parameter e and the Bogoliubov parameter \/3\. In terms of the reduced variables x 2 , 
X3 we have 

pmodin/. „ „ x 1 f 1 ~ COs[fci?7o(l + X 2 - X 3 )\ 1 - CQ S [fcl?7o(l - X 2 + X 3 )] , 

b {kir] ,x 2 ,x 3 ) = < r 1 r h 

X 2 X3 { X$(l + X 2 - X3) X%(1 - X 2 + X3) 



1 - cos[/ci?7o(-l + x 2 + x 3 )] 
(-1 + x 2 + x 3 ) 



(3.3) 



We now explicitly see the dependence of the shape function on |fei7/o| = ^^°D = M/H, 
namely the ratio of the physical cut-off scale to the Hubble parameter, as was explained 
earlier. Again, not fixing the combination |/cir/o| to be /ci-independent results in an obvious 
breaking of scale-invariance and would not allow us to use the introduced tools for compar- 
ison with the available templates. The cut-off scale should be significantly larger than the 
Hubble scale and we will typically consider it to be somewhere in between 10 2 — 10 3 . This 
implies the shape function is rapidly oscillating, which complicates the evaluation of the 
integrals to determine the cosines and fudge factors with the available templates. When 
possible the integrals were evaluated analytically in the limit |A;i?7o| 3> 1. Let us compute 
the squared norm of the modified initial-state shape function given by 



i ri 



2 



4 4 
X 2 X 3 



F moam (k lV0 ,x 2 ,x 3 ) = dx 2 dx 3 F™ am (k m ,x 2 ,x 3 ) 

JO Jl-x 1 J 

= \hr)o\ + ^loglfci^ol + 6.05 , (3.4) 
bU 4 

where the integral in x 2 , x 3 space is over the (triangle) domain < x 2 < 1, 1 — x 2 < x 3 < 1. 
From Eq. (3.4) we can derive some important conclusions about the detectability of this 
non-Gaussian signal. In an ideal situation the data analysis would be performed using 
the theoretical template Eq. (3.2) to directly infer on the non-Gaussian amplitude. As 
previously discussed this is the product of the normalization A times the norm of the 
shape function, A\F modm y Hence, the best one can do by using an observational template 
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perfectly aligned with the theoretical prediction is a leading enhancement factor of order 
v/|fci?7o[. However, such an enhancement is lost when the data analysis is performed using 
a local template, due to the integrated nature of the non-Gaussian analysis. Evaluating 
the scalar product as defined by Eq. (2.4) between the initial-state modification template 
and the local one we find 

F™ din (k lVo ,x 2 ,x 3 ) ■ F local (x 2 ,x 3 ) = 
f dx 2 [ dx 3 F modm (k m ,x 2 ,x 3 )F loca \x 2 ,x 3 )x 4 2 x 4 3 = 3log\k im \ + 18.96. (3.5) 

JO Jl-x 

We can already conclude that when using the local template to probe modified initial state 
non-Gaussianities the enhancement factor is further reduced to become only logarithmic 
in |fci?7o|, instead of the -v/|fci?7o| enhancement that can be achieved in the optimal case. It 
is worth remarking that the constant parts in the results, for both the norm and the scalar 
product, depend on how some singular integrals, those independent of kirjo, are being cut- 
off. The singular integrals always blow up in the local (squeezed) limit, corresponding to 
one of the momenta being much smaller than the other two. Fortunately, a natural cut-off 
is given by the fact that only a finite range of modes contribute to the CMB. Specifically, 
the ratio between the smallest and the largest observable scales on the CMB is roughly 
equal to 10 -3 , which is used to regularize the local type integrals. Throughout this paper 
we will always quote results using this cut-off if necessary. For the squared norm of the 
local shape, which clearly blows up in the squeezed limit, the need for this cut-off is most 
apparent. To be more explicit, using the 10~ 3 cut-off the result for the squared norm of 
the local shape equals |i ?local (x2, x 3 ) | = 176.5. From the scalar product and the local and 
modified initial-state norms we can infer the cosine factor, which reads as 



Cos 



F modin Recall = F™**(k m , X 2 , X 3 ) • F^\x 2 ,X 3 ) 

\ \F* aodi *(k 1 rio,X2,x 3 ) \ \F loca \x 2 ,x 3 )\ 



7.53 • 10- , (18-96 + 3bg 1*^1) (3 6) 

6.05 + ^\kii]o\ + f log \kirj \ 

The fudge factor necessary to transform the limits on local type non-Gaussianities into 
constraints on modified initial-state non-Gaussianities is given by 

Ap T^modin Recall = F^(k lVo , X 2 , X 3 ) ; F^[x 2 , X 3 ) 

= 5.67 -10~ 3 (18.96 + 3 log \kii ]0 \) . (3.7) 

As we already concluded from the scalar product alone, the fudge factor is logarithmically 
dependent on kirjQ. In addition, the coefficient is also relatively small, implying that over a 
realistic range range of \kir]o\ values, the fudge factor can essentially be considered constant. 

Similarly, we calculate the scalar product between the equilateral template and the 
modified initial-state distribution, which a priori can be expected to depend on k\r]Q as 
well. Surprisingly, the leading k\rjo dependent terms cancel and the only contribution 
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Figure 1: Cosines factors between the initial-state modification shape and the local (solid line), 
equilateral (long dashed line) and the enfolded template proposal (short dashed line) as functions 
of \kir] \ in 3-D. 

comes from a, cut-off independent, constant number for the scalar product 

F™ di *(k lVo ,x 2 ,x 3 ) ■ F e ^ il (x 2 ,x 3 ) = 
dx 2 [ dx 3 F modm (k lVo ,x 2 ,x 3 )F C(inil (x 2 ,x 3 )xixj = 6.5. (3.8) 

Jl-x 

We can therefore conclude that all enhancement is lost when using the equilateral template 
to probe modified initial-state non-Gaussianities. As for the local template this will imply 
a constant fudge factor, even though the theoretical non-Gaussian distribution is linearly 
enhanced in enfolded triangles. The squared norm of the equilateral shape functions is also 
cut-off independent (i.e. finite), and the numerical integration gives \F e ^ uil (x 2 ,x 3 )\ = 7.9. 
Combining this with the squared norm of the modified initial-state shape function this 
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leads to the following expression for the normalization independent cosine 



Cos 



^imodin ^iequil 



F modhl (k lV0 , X2, x 3 ) • F e ^(x 2 , x 3 ) 
| F modin( /cir/0) X2j X3 )| |i?equil( X25 

2.31 , . 

^^^=^^^= . (3.9) 

6.05 + gg|/ci7? | + | log \klT]o\ 

Even though the cosine is a function of kirjo, while the scalar product is not, the fudge 
factor will also be independent and equals 

F mod ^(k m ,x 2 ,x 3 ) ■ F e ^(x 2 ,x 3 ) 



A, 



^imodin ^equil 



\F^(x 2 ,x 3 )\ 2 
0.82. (3.10) 



The constancy of the fudge factor explicitly confirms that all enhancement due to the large 
\kiT]o | parameter is lost. In Figure 1 and 2 we plot the cosine and fudge factors between 
the initial-state modification and the local (solid line) and equilateral (long dashed line) 
templates as function of k\T]o. From the plot of the cosine factor we see that indeed the 
local and equilateral templates poorly overlap with the modified initial-state distribution 
as the cos [F modin ,F x ] < 0.6. We conclude that although the non-Gaussian amplitude of 
initial-state modifications is linearly enhanced in enfolded triangles, the measured local and 
equilateral templates are completely insensitive to this localized enhancement, thus spoiling 
any chance of obtaining a stringent bound on departures from the standard Bunch-Davies 
vacuum state. 

Consequently, probing standard slow-roll modified initial-state non-Gaussianities is 
impossible unless a new template distribution is introduced which, unlike the local and 
equilateral templates, is sensitive to the localized enfolded enhancement. As pointed out 
in the previous discussion, using a perfect template will lead to a signal enhancement of 
y[fei77o[. In Section 5 we will describe a first proposal for such an improved, factorized, 
enfolded template. In the next section we will focus on the combined effect of a specific 
higher derivative correction and an initial-state modification on the bispectrum. 



4. Adding higher derivative corrections 

As was shown in the previous section, the enhancement effect of an initial-state modification 
in the bispectrum, assuming standard slow-roll inflation, is impossible to probe using the 
currently available local or equilateral templates. What we would like to study is whether 
the same conclusion holds after adding higher derivative corrections, which according to [54] 
could be even more sensitive to initial-state modifications. A priori, one might expect 
similar conclusions, that even though there is a strong enhancement effect in the enfolded 
triangle limit, its measure in the space of all triangles versus the local or equilateral template 
will again be too small to allow detection. We consider the addition of a dimension 8 higher 
derivative term to the scalar field lagrangian of the following form: 
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Figure 2: Fudge factors as functions of |fci?7o| for the templates as in figure 1. 

where the scale M corresponds to the high energy cut-off scale and 'natural' corrections 
would correspond to a coupling A ~ 1. This higher derivative correction is the same as the 
one discussed in [54] and was first studied in [19]. We provide a detailed derivation of the 
corresponding bispectrum in Appendix B. As shown in [19], assuming the standard Bunch- 
Davies vacuum, it leads to non-Gaussianities of the equilateral type with an amplitude 
f^ 11 oc ( ^m^ ^ ^ e which, at best, can be of order 1 (in order not to spoil the higher 
derivative expansion) . In the interaction Hamiltonian for the relevant perturbation variable 
£ it leads to an additional term of the form 

AH> = ~ j d* X a(r,) U) C' (C 2 - (d lC f) . (4.2) 

As first shown in [54], and repeated here in the appendix, the associated bispectrum cor- 
rection due to an initial-state modification is a complicated function of the comoving mo- 
menta. Most importantly, compared to the result obtained for the standard slow-roll 
computation, after integrating over conformal time one now finds terms proportional to 
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r/o and t/q, in addition to contributions independent of r/o- The different powers of 7/0 can, 
as before, be combined with one of the comoving momenta k\ to give the large number 
|&l?7o| = M/H 3> 1. Consequently, the amplitude of the three-point function is expected 
to be dominated by contributions proportional to t/q. Collecting the leading contributions 
and neglecting terms that are not (locally) enhanced at t/q order, we obtain 



where 



(CfciCfc 2 0c3)nBD 



1 - cos(fc i r ?0 ) sin(fcj7/ ) , 

7]q ft | V{kj,k j+ i,kj +2 ) 



-7/0 (cos(%r/ ) Q(kj,k j+1 , k j+2 ) 



(4.3) 
(4.4) 



V(k j ,k j+1 ,k j+2 ) 



Q(kj,k j+1 ,k j+2 ) 



-4k j+ ik j+2 (k j+1 + k j+2 )(k 2 j+1 + fc| +2 + k 



•j+ikj +2 ) 



+2 x Yl H k t ~ 4k j+1 k j+2 ). 

i 



In the above expression kt = Yli k i represents the sum of all (absolute values of) comoving 
momenta and kj is defined as before kj = kt — 2k j, i.e. the sum of comoving momenta with 
one of the signs reversed. The terms in Eq. (4.3) are enhanced to an order r/g only in the 
collinear limit, while being suppressed for all other triangular configurations. In contrast 
the term in Eq. (4.4) is enhanced by t/q over the full triangle domain and is therefore 
expected to be the dominant contribution. This was apparently not noticed in [54] , maybe 
because the collinear limit was assumed from the start. As a result the shape of the higher 
derivative bispectrum with a modified initial-state is not of the expected enfolded type. 
To be explicit let us rewrite the dominant overall enhanced contribution in terms of the 
two variables x 2 = jj^ and X3 = where we scaled out the usual A;^ 6 and absorbed the 
enhancement factor (fcir/o) 2 into the non-Gaussian amplitude /nl> producing the following 
shape function 



F 



HD-dom 



(kirj ,x 2 ,x 3 ) 



1 



{cos(fci?7o(x2 + x 3 - 1)) [(1 



x 2 + x 3 J 



4x 2 x 3 



+ cos(kii] (-x 2 + x 3 + 1)) [(1 + x 2 + x 3 ) 2 - 4x 3 ]) 

+ cos(A;i 7/0(2:2 - x 3 + 1)) [(1 + x 2 + x 3 ) 2 - Ax 2 ] } . (4.5) 

As before, scale-invariance of the bispectrum therefore requires the combination |fci7/o| = 
M/H ^ 1 to be k\ independent. Note that the cosines appearing in this shape function 
imply that the three-point function is constantly changing sign. The norm of the full 
higher derivative modified initial-state distribution is well approximated using only the 
contribution described by the shape function ^ HD - dom . We find, after averaging over the 
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cosine, that |i?HDmodin|2 ^ j^HD-domp _ 23.34 Since it is the normalization independent 
combination /nl \ F\ that is actually being measured, we conclude that a perfect template 
would be sensitive to the full {k\rjo) 2 enhancement factor. Because this number could 
be as large as 10 6 it indicates that the higher derivative terms are extremely sensitive to 
initial-state modifications, potentially leading to strong constraints on departures from the 
standard Bunch-Davies vacuum state. This derives from the fact that the nonlinear higher 
derivative interaction plays an important role at sub-horizon scales. Sub-horizon particle 
occupation numbers as a consequence of the modified initial-state allow for the generation 
of a significant non-Gaussian signal due to the crucial presence of the higher derivative 
interactions at that stage. This is different from the standard slow-roll situation where the 
required (gravitational) nonlinearities are far less important at sub-horizon scales. 

Unfortunately though, due to the oscillating sign nature of the dominant contribution 
Eq. (4.4) the currently available local and equilateral templates are extremely insensitive 
to this term, i.e. the scalar products between ^? HD - dom anc l the equilateral and local 
templates are suppressed because of cancellations inside the scalar product integral. A 
quick inspection of the scalar product integral reveals it could scale as l/(/cir/o) times some 
oscillating function of k\rjo, which would reduce the overall (kirjo) 2 level of enhancement 
by at least one power. It is for this reason that we have kept the locally enhanced terms 
of Eq. (4.3), since these could give rise to contributions in the scalar product of similar 
order in kirjQ. The sine term in Eq. (4.3) is overall enhanced with one power of k\rjo and, 
based on the results in the previous section, the localized enhancement due to the single 
kj in the denominator is expected to disappear after calculating the scalar product with 
the local or equilateral template, neglecting possible logarithmic terms. The cosine term is 
locally enhanced by a factor (kirjo) 2 due to the squared kj dependence in the denominator. 
One of those enhancement factors is again expected to be lost after performing the scalar 
product integral with the local or equilateral templates. 

We anticipate a linear kirjo scaling at best (neglecting possible logarithmic terms) and 
we should keep track of all terms in Eq. (4.3) and Eq. (4.4) when evaluating the scalar 
product with the local or equilateral template. Unlike the previous section we were unable 
to perform the relevant integrals analytically and instead relied on a numerical approach, 
fitting the scalar product integral results for a large sample of kirjo values to estimate the 
k±rjo dependence. The relevant shape function is identified in exactly the same way as the 
dominant contribution i? HD - dom 5 except this time no overall factors of kirjo are absorbed 
into the definition of the non-Gaussian amplitude 

As before we assume that the Bogoliubov parameter lZe(f3kj) ~ Having fixed the non- 
Gaussian amplitude and shape function ambiguities we find that the leading order behavior 
of the scalar product with the local template distribution is well approximated by 

F local • F HDmodin « (ki V0 ) (-72 + lOlog \hrj \) , (4.6) 

4 Remember that the overall (kiijo) 2 enhancement factor was absorbed into the non-Gaussian amplitude 
/nl, explaining why it does not show up in the norm. 
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confirming the general expectation on the order of magnitude of the result. We should 
point out that the relative minus sign between the different terms in Eq. (4.3) is the source 
of the relative minus sign in the final result Eq. (4.6). Rather unfortunately, the different 
coefficients conspire in such a way that the scalar product has a minimum and then crosses 
through zero in the domain of interest 100 < |&i?7o| < 1000. This implies significantly 
smaller fudge factors, for a small range of |&i??o| values, than would be expected on the 
basis of scaling alone. Another consequence of this, confirmed by the numerical results, is 
that the unavoidable oscillatory contributions, that we neglected when fitting the numerical 
data to produce Eq. (4.6), are bound to give rise to relatively large corrections in the \kirjo\ 
domain of interest. The makes the full structure of the scalar product rather complicated. 
Although the general trend is nicely described by a linear plus logarithmic scaling with 
|/ci?7o|, in the |fci?7o| domain of interest the actual value of the scalar product fluctuates and 
can deviate from the expected order of magnitude for some values of |A;ir/o|. Dividing the 
scalar product by the norm of the local distribution (which was already computed in the 
previous section) one obtains the fudge factor 

\F^(x 2 ,x 3 )f 
5.7 • lO" 3 ^!^! (—72 + lOlog |fci7? |) • (4.7) 

We conclude that a |fci?7o| log \ kiVo I enhancement remains, which is almost one power of k\r]o 
less as compared to the optimal scenario. For the fudge factor, as for the scalar product, 
the same cautionary remarks apply. The above result describes the average trend and the 
detailed numerical results show that fluctuations can have a significant effect on the actual 
value of the fudge factor in the |fci?7o| domain of interest. As we will see, the results for the 
2d fudge factor with the local template exhibit a similar complicated behavior as a function 
of kirjo, although the actual numbers for the fudge factor, due to the larger coefficients, are 
roughly one order of magnitude larger. For the equilateral template the final scaling result 
is the same, although somewhat surprisingly, it is the dominant contribution in Eq. (4.4) 
that is solely responsible for the final result. As one can easily check, both terms in Eq. (4.3) 
are in fact maximizing exactly at the line x% + £3 — 1 = 0, whereas the equilateral template 
is exactly vanishing at the line X2 + 23 — 1 = 0. The result is a suppressed contribution to 
the scalar product which is negligible in the limit of large | A^i 770 1 as compared to the other 
contribution. Incidentally this observable enhancement is the same as reported in [54], 
but the underlying reason is very different. It is a consequence of using a template that 
is far from optimal and it should be possible to achieve significantly higher sensitivity by 
constructing a more suitable template to analyze the data. In particular, the non-Gaussian 
signal described by Eq. (4.3) and Eq. (4.4) is not of the enfolded type and has strong 
oscillatory features, which might allow for a clear distinction from other non-Gaussian 
sources. 

The generic appearance of at least a single factor of |fci??o| = M/H in the fudge fac- 
tor with respect to the local (or equilateral) template implies an enhancement possibly as 
large as 10 3 , ignoring the fluctuations of the fudge factor as a function of /ci?7o- At the 
start of this section we mentioned that in the standard Bunch-Davies vacuum the higher 
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^iHDmodin ^local 
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derivative term would give rise to a maximal of order 1. Compared to the original 

higher derivative non-Gaussian amplitude, the modified initial-state amplitude Eq. (4.6) 
introduces an additional suppression with the Bogoliubov parameter. On the other hand 
the fudge factor introduces a linear |/ciryo| = M/H factor enhancing the original higher 
derivative non-Gaussian amplitude by (3 (M/H) when probed with the local or equilateral 
template. The CMB two-point power spectrum constrains the Bogoliubov parameter al- 
ready at the 10 -2 level, so at best this would allow for a local or equilateral non-Gaussian 
amplitude of order 10 due to initial-state modifications, assuming M/H ~ 10 3 . This might 
be detectable in the future, although there are many other sources for a local or equi- 
lateral non-Gaussian signal at that level. In section 6 we will confirm the same level of 
enhancement by computing the projected 2d fudge factor and use the most recent WMAP 
constraints on local type non-Gaussianities to derive an order of magnitude constraint on 
the Bogoliubov parameter. 

5. An enfolded template proposal 




Figure 3: From left to right: a squeezed, equilateral and squashed triangle. 

In the absence of higher derivative corrections we have shown that a non-Gaussian 
signal due to a modified initial-state, which maximizes in collinear triangles, cannot be 
probed using the available local and equilateral templates. Both templates are not sensi- 
tive enough to the localized enfolded enhancement to give rise to a significant (preferably 
power law) dependence of the fudge factor for large |&i?7o| = M/H. Instead the local and 
equilateral fudge factors are at best logarithmically dependent on |fcir/o|. To see if one can 
improve on that situation one would like to introduce a more suitable template, one better 
aligned with the theoretical prediction of modified initial-state non-Gaussianities. 

The distinguishing feature of such a template should be that it maximizes in the en- 
folded (or squashed) triangle limit, as opposed to the local and equilateral templates which 
maximize in squeezed and equilateral triangles respectively. Whereas the squeezed trian- 
gle is obtained by taking one of the comoving momenta to zero ki — ► 0, and equilateral 



-17- 



corresponds to all momenta equal, enfolded triangles imply two collinear momenta, and 
therefore ki = kj + k m with i ^ j ^ m and kj and k m representing the two collinear 
momenta. Clearly the squeezed, equilateral and squashed triangle limits exhaust all pos- 
sibilities, which are shown in figure 3, and nicely correspond to the three different classes 
of theoretical non-Gaussian predictions: local, equilateral and enfolded. So besides its 
potential theoretical relevance, also from the point of view of completeness it might be 
worthwhile to develop a third factorized template shape that would maximize in the en- 
folded triangle limit. This would introduce a third non-Gaussian observable /j-gj measuring 
the enfolded amplitude. 

The reason why one cannot directly compare theoretical predictions to the CMB data 
and needs especially designed templates lies in the computational complexity of recon- 
structing the non-Gaussian amplitude from the two-dimensional CMB temperature data. 
The projection of a 3-point correlator to a 3-point function in spherical harmonic space 
involves the Wigner 3j symbol (to construct the angular averaged bispectrum) and a com- 
plicated integral over transfer and Bessel functions. This is computationally very challeng- 
ing, scaling as iV 5 / 2 ; N 1 / 2 for every multipole / and N for the averaging over m, where 
N equals the total number of pixels in the CMB map. In the last few years different sug- 
gestions have been made to accomplish a reduction of computational time [33,39-46]. A 
significant reduction in the number of calculations can be achieved if the three-point func- 
tion is factorizable in its momentum dependence, i.e. F(k\, ki, ^3) — ► f a (ki)fb(k2)fc(k3), 
leading to a reduction from iV 5 / 2 to iV 3 / 2 . As it turns out, local type non-Gaussianities 
are indeed described by a factorized shape function -F(fci, &2, £3), whereas the theoretical 
predictions for non-Gaussianities of the equilateral and enfolded type are not factorizable. 
This makes the direct comparison of equilateral and enfolded type non-Gaussianities to the 
two-dimensional CMB data extremely difficult for the time being, although recently some 
progress has been made to allow for a more direct comparison of arbitrary signals [56-58] . 
Up to now one instead relies on factorized template approximations to the theoretical 
signals. Not so long ago an equilateral template has been successfully identified and com- 
pared to the data [34,35,37,43], but an enfolded observational template has not yet been 
constructed. Below we will construct a first proposal for a factorized enfolded template 
and analyze how much better it will be able to constrain the modified initial-state non- 
Gaussianities discussed in section 3. 

Looking at Eq. (3.2) it is clear that the three-point correlation function due to initial- 
state modifications is not factorizable. As explained this non-Gaussian shape function is 
the result of adding a minus sign to one of the comoving momenta to first order in the 
Bogoliubov parameter /3k [14,54]. For instance a term behaving as l/(fci + k<i + ks) would 
change to 1/ (— k\ + k% + k^) plus permutations. This suggests an approach where one starts 
with the factorized equilateral template eq.(2.9) and just replaces ki — > —ki, symmetriz- 
ing over all the indices. Applying this idea produces the shape function F{k\, k2, k^) = 
—F equ i\(ki, &2, ^3), which does not yet resemble the desired enfolded distribution nor does it 
add additional information, i.e. it is simply the equilateral shape multiplied by a minus sign. 
Fortunately though it requires only a small modification to come up with a factorized shape 
function that seems to be a reasonable candidate for an enfolded template. Starting from 
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Figure 4: The enfolded template shape F(x2,xs) 



the equilateral shape function, replacing k — > —k, introducing X2 = ky/kx and X3 = k$/k\ 
and plotting the obtained distribution F(x2,xz) times the appropriate measure factor 5 



x 2 x 3 



, it becomes apparent that a term proportional to l/kfk^k^ acts as a kind of constant 
'normalization' of the template. Additional \ jk\k\k\ terms therefore simply lift or lower 
the whole graph. By adjusting the number of such terms, so F = —F eqml + cjk\k\k\, 
it is possible to lift the obtained shape in such a way that it resembles an enfolded type 
distribution, maximizing on the line ky, + k$ ~ fei, corresponding to enfolded triangles. We 
find that the best choice requires adding only one such term to the — -F e quii distribution, i.e., 
c = 1 (see Appendix A). Consequently our proposal for the factorized enfolded template, 
as a function of x%, X3, becomes 



F enf (x 2 ,x 3 ) = 6 



x 



1 3/1 
-3+2 perm H 2 — ( 2 3 + " P erm 



(5.1) 



We have plotted the template shape function in figure 4. In appendix A we explain in 
what sense c = 1 corresponds to the optimal choice and the details are presented on 
how to translate this template into an observable using the 'fast best estimator' approach 
developed in [33,39,43,44,46]. 

We should determine how well the proposed enfolded template overlaps with the the- 
oretical modified initial-state three-point function of Eq. (3.2). To quantify this we will 



3 This is the relevant quantity because of the measure in the scalar product of eq. (2.4) 
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perform the same analysis as in section 3, calculating the scalar product, cosine and fudge 
factor, now using the enfolded template. The squared norm of the modified initial-state 
shape function is given by Eq. (3.4) and a numerical integration gives \F eni (x2, x 3 )\ 2 = 4.34. 
After computing the scalar product, this leads to the following expression for the cosine as 
a function of | Ari 770 1 



Cos 



^modin ^cnf 



F™ d ™(k im ,x 2 ,x 3 ) ■ F™ f (x 2 ,x 3 ) 
I F modin (A;i 7io,x 2 ,x 3 )\ I F enf (x 2 ,x 3 )\ 
(2.80 + log \ki Vo \ 



4.80 -IP' 1 — v 01 WU|/ (5.2) 



6.05 + ^Mol + § log 



We have plotted this function in figure 1, including the cosine functions between the mod- 
ified initial-state shape distribution and the local and equilateral templates. As should be 
clear from the plot the cosine between the enfolded template and the theoretical distribu- 
tion is closer to one, but there is certainly room for improvement. As the parameter |fci?/o| 
grows, the enfolded template will depart more from the theoretically predicted modified 
initial-state shape. Nevertheless, the enfolded template has significantly higher overlap 
with the theoretical distribution than the local or equilateral templates, at least in the 
comoving (3-D) momentum space. Since the cosine is the relevant quantity that deter- 
mines the relative improvement, comparing to the plots (see also table 1) for the local 
and equilateral cosines one concludes that a rough 35 — 45 percent level of improvement 
should theoretically be achievable using the enfolded template. This is certainly not enough 
to derive interesting constraints for the theoretically predicted modified initial-state non- 
Gaussianities, as can be seen more directly by turning our attention to the fudge factor, 
which explicitly identifies the level of kiT]o enhancement. From the expression of the cosine 
it is straightforward to read off the fudge factor 



A; 



^imodin jpcnf 



0.65 + 0.23 log \kmo\ (5.3) 



which disappointingly implies that the enhancement, or the sensitivity, that would be 
achieved using the proposed enfolded template is still only logarithmic in | /c 1 770 1 , far removed 
from the maximally attainable level of y[feir?o| enhancement. Comparing to the local and 
equilateral fudge factors the coefficients are bigger which means the enfolded fudge factor 
will be larger and changes considerably over the natural range of | fci 770 1 (10 2 — 10 3 ). This is 
however still far removed from the (power law) level of enhancement that one would need 
to derive interesting constraints on the Bogoliubov parameter from the bispectrum data. 

In tables 1 and 2 we have collected the cosine and fudge factor results for the dif- 
ferent templates with respect to each other and some important theoretical predictions 
(equilateral higher derivative and enfolded modified initial-state). Note that the modified 
initial-state entries typically depend on the (large) parameter kirjo, which for the cosines 
has been denoted by the range of possible values, whereas for the fudge factors we have 
chosen to explicitly write the function. Independent from the original theoretical moti- 
vation, one could argue that the enfolded template nicely completes a general analysis 
of non-Gaussian signals. Table 1 then shows how much complementary information each 
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Local 


1 


0.41 


0.68 


Equil. 


0.41 


1 


0.49 


Enf. 


0.68 


0.49 


1 


HD 


0.45 


0.99 


0.59 


Mod 


0.6 - 0.3 


0.6 - 0.4 


0.9-0.6 



Table 1: The 3d Cosine 



Shape Fy 


3-D Fudge (F x 


= F\ oc ) 


3-D Fudge (F x = F eq .) 


3-D Fudge (F x 


= F enf ) 


Local 


1 




1.94 


4.29 




Equil. 


0.09 




1 


0.66 




Enf. 


0.11 




0.36 


1 




HD 


0.10 




1.07 


0.86 




Mod 


0.11 + 0.0171o£ 


? \hr]o\ 


0.82 


0.65 + 0.23 log 


\hr]o\ 



Table 2: The 3d Fudge Factors. Note that HD are the higher derivative contributions from [40]. 
We added these for completeness and to show consistency with the results in [40] . 



template would provide. Compare this to a decomposition of a general vector into a set of 
basis vectors. Ideally, one would prefer to come up with a set of orthogonal basis shapes. 
Instead, the local, equilateral and enfolded template are far from orthogonal, but each 
does provide complimentary information that can be precisely quantified in terms of the 
different cosine values 6 listed in the table. Decomposing a general three-point signal in 
these template shapes might therefore still be useful, even though the enfolded template 
by itself is unable to probe modified initial-state non-Gaussianities. 

6. Two-dimensional bispectrum results 

Differently from the 3-D case the CMB temperature anisotropies are a 2-D projection of the 
linearly evolved primordial curvature perturbation field, hence the 2-D bispectrum is the 
result of the convolution of the shape function with the photon transfer function projected 
on the sky. In the following we give a brief review of the basic formalism, then we will 
discuss the results of the CMB bispectrum computation. 

Let us consider the standard spherical harmonic decomposition of the CMB tempera- 
ture fluctuation along the direction h of the sky, 

— (n) = 5>riT(n), (6.1) 

l,m 

the multipole coefficients af 1 contain all statistical information about the temperature 
anisotropy field, and are the starting point to construct the various correlator functions. 

6 Note that since we are dealing with a function space, obviously a complete decomposition would formally 
require an infinite set of basis functions on the relevant triangle domain. 
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The angular bispectrum in multipole space is given by 

KLTr 3 =< <CCC >. ( 6 - 2 ) 

and assuming rotational invariance, the angle-averaged bispectrum reads as 

*•<*- E [!!!)<<'«>■ <«*) 



mi,m2,rrt3 



mi m2 ms 



Substituting the expression of the multipole coefficients in terms of the photon transfer 
function and the primordial curvature perturbation, Eq. (6.3) becomes 



with 6^/3/3 the reduced bispectrum given by, 



2 X3 



7T 



b hhh = [ - ) I dxdk 1 dk 2 dk 3 (xkik2k3) 2 ji 1 (k 1 x)ji 2 (k2x)ji 3 (k 3 x 



xF(ki,k 2 ,k 3 )A ll (k 1 )A h (k 2 )A l3 (ks), (6.5) 



where A;. {ki) is the photon transfer function and ji^kix) is the Bessel function (for a 
detailed derivation see [22]). As discussed in [40], the evaluation of the reduced bispec- 
trum is computationally challenging, on the other hand in the flat-sky approximation the 
computation is simplified, since for example the integral over the Bessel function does not 
explicitly appear. In such a case one has: 

Khh = (r ° (2 ~^ )2 J dkfdkldk z 3 6(kt + kl + k z 3 )F(k[,k' 2 ,k' 3 )^ 

(6.6) 

where to and tr are the conformal time today and at decoupling respectively, kf is the 
component of the wave-vector orthogonal to the plane tangent to the last scattering sur- 
face, and k! = yj (k z ) 2 + I 2 / (tq — tr) 2 ; the photon transfer function along the orthogonal 
direction is given by 

Mk z ) = J j^y 2 S{^m 2 + I 2 /(r - T) 2 y k *\ (6.7) 

where S{...) is the CMB source function. The presence of the delta-function ensures that the 
projected modes form closed triangles. Hereafter we assume a vanilla LCDM cosmology 
with model parameter values corresponding to the WMAP 5-years best-fit model. All 
cosmologically relevant quantities such as the source functions have been computed with the 
publicly available CMBFAST code [59]. We then have evaluated the reduced bispectrum 
for different shapes using Eq. (6.6). 

Following [40] we introduce a scalar product 

Bx-B Y = f B * l2 fl h r . ( 6 - 8 ) 
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where fi x i 2 i 3 is a combinatorial factor which is 1 if all three multipoles are different, 2 if two 
of them are equal and 6 if all of them are equal; C\ is the angular CMB power spectrum 
which includes the experimental noise evaluated in the Gaussian approximation assuming 
WMAP experimental characteristics. The cosine reads as 

Cos™ [B X ,B Y ] = (6-9) 

and the fudge factor as 

A 2 F D [Bx,B Y )= * x B ^ . (6.10) 

"local - "local 

In principle we may expect that the non-Gaussian signal given by different triangular 
shapes on the CMB differs from that obtained in 3-D, since triangles of different shapes in 3- 
D can be projected into the same 2-D configuration. If this is the case then the values of the 
cosine factors should be shifted upwards. Evaluating the cosine and fudge factors between 
local and equilateral shapes we find Cos 2D [F local , F c ^ 1 } = 0.62 and A|P[F local , F equil ] = 
0.13 respectively, which is consistent with the results presented in [40]. The cosine and 
fudge factors between the initial-state modification shape function and the local, equilateral 
and enfolded templates are shown in figure 5 and 6 as a function of the kirjo parameter 
respectively. 

We may notice a trend similar to that inferred from the 3-D evaluation. In particular, 
the cosine decreases as a function of kirjo for all three templates, whereas the fudge factors 
are constant for the local and equilateral case, and increasing for the enfolded template. 
The enfolded template has the largest overlap with the initial state modification shape, 
although not significantly better then the local one. Overall the cosine values are slightly 
larger than what we have found in the 3-D calculation. This is because different trian- 
gular configuration in 3-D can be degenerate in the 2-D, hence the projection tends to 
systematically increase the overlapping between different templates. 

We have also computed the reduced CMB bispectrum for the initial-state modification 
in the presence of higher order derivative terms as given by the shape function Eq. (4.3) and 
Eq. (4.4). We confirm the enhancement effect discussed in section 4, as an example calcu- 
lating the fudge factor with the local template for kirjo = 100 we find |A| D | ps 100, wher cas 
for kirjo = 10 3 we obtain |A|P| ~ 6000. This clearly shows that the non-Gaussianity in- 
duced by initial state modifications is enhanced by the presence of higher derivative terms 
leading to potentially large detectable signals. The specific functional dependence of the 
fudge factor on kir/Q in the range of interests (10 2 — 10 3 ) is far from trivial due to the 
interplay of the different terms in Eq. (4.3) and Eq. (4.4). These contain oscillatory factors 
that leads to a modulated oscillations of the fudge factor dependence on kirjo. Besides, 
the same term causes an oscillatory dependence of the reduced bispectrum as function of 
the multipoles. These oscillations are responsible for cancellations in the evaluation of the 
cosine factor, hence leading to a very small overlap with the other (non-oscillatory) tem- 
plates. For example evaluating the cosine factor with the local shape for kirjo in the range 
10 2 — 10 3 we find |CosfP| ~ 0.01, implying that current observational templates are not apt 
to detect such a non-Gaussian signal. The determination of an optimal observational tem- 
plate that can account for the feature produced by this type of non-Gaussianity is therefore 
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Figure 5: 2-D cosine factors between the initial state modification template and the local (solid 
line), equilateral (long dashed line) and enfolded (short dashed line). 

necessary and we leave this search to future work. Nevertheless using the above estimates 
of the fudge factor, we can use current constraints on local non-Gaussianity to infer limits 
on /NL >m ° dm an d thus on the Bogoliubov parameter. As an example the prediction for a 
local non-Gaussian contribution due to a modified initial-state in the presence of a higher 
derivative operator reads 



|A/i 



local I 
NL I 



6 



M 2 Pl H 2 
M 4 



IA 



2D I 



(6-11) 



where we assumed that the coefficient A in Eq. (4.6) equals one. Using the observed am- 
plitude of the power spectrum, the slow-roll parameter e can be replaced with jjfir- 
Assuming M/H ~ 10 3 this gives rise to the following constraint on the Bogoliubov pa- 
rameter, using the latest WMAP 5-years upper limit on an d the result for the fudge 
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Figure 6: 2-D fudge factors between the initial state modification template and the local (solid 
line), equilateral (long dashed line) and enfolded (short dashed line). 



factor \Af\ « 6000 



4 • 10 3 \p\ < 111 



(6.12) 



Surprisingly, this corresponds to a relatively strong bound on the Bogoliubov parameter 
\P\ < 3 • 10~ 2 . Note that in most proposals one expects the Bogoliubov parameter to be 
a function of H/M. For instance in the New Physics Hypersurface scenario \f3\ is predicted 
to be linear in H/M [27, 28] . Using the results above for H/M ~ 10 -3 this predicts at 
best an order 10 contribution to the local non-Gaussian signal, which could increase to an 
order 10 2 contribution for H/M ~ 10~ 2 . We should stress that an ideal template could 
improve the limits on the Bogoliubov parameter by another factor of M/H. 
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7. Conclusion 



We have analyzed inflationary three-point correlators as a result of a small departure from 
the standard Bunch-Davies vacuum. In the simplest scenario where we avoided higher 
derivative interactions we confirmed that the initial-state modification causes the three- 
point correlator to maximize in the collinear or enfolded triangle limit, corresponding to a 
uniquely different shape as compared to local and equilateral non-Gaussian signals. Since 
the maximal signal scales linearly with the cut-off scale M, the non-Gaussian amplitude 
in enfolded triangles can be quite large, perhaps allowing for detection or providing in- 
teresting constraints on departures from the Bunch-Davies vacuum. However, by com- 
puting the scalar products, and consequently the cosine and fudge factors, between the 
theoretical prediction for the three-point function and the existing local and equilateral 
observational templates, we concluded that essentially all enhancement is lost due to the 
inefficiency of the available observational templates combined with the localized nature of 
the enhancement. Although [54] reached a similar conclusion, their argument was very 
different, relying on the projection to the two-dimensional CMB sphere. Instead, we have 
shown that the currently available method of comparing theoretical three-point functions 
to CMB bispectra, involving observational templates and the necessary integration over all 
triangles, already removes most sensitivity to localized enhancements in enfolded triangles, 
even before projecting to the two-dimensional CMB sphere. 

The situation can in principle be improved by constructing a suitable enfolded tem- 
plate (better) adapted to the theoretical prediction. Moreover, from a general non-Gaussian 
analysis point of view, the introduction of an enfolded template might be interesting in 
itself, potentially providing complementary information in addition to the local and equi- 
lateral templates. The enfolded template proposed here was unfortunately only a marginal 
improvement over the local and equilateral template, still being insensitive to the local- 
ized enhancement. It would certainly be worthwhile to look for a more optimal enfolded 
template that could approach the theoretically maximum level of sensitivity to the local- 
ized enhancement, corresponding to a = \J M/H dependence of the corresponding 
fudge factor. 

After adding a specific higher derivative term we surprisingly found that the localized 
nature of the enhancement is substituted by an overall enhancement of (/ci%) 2 that can 
be absorbed directly into the non-Gaussian amplitude /nl- Even though (sub-leading) 
terms exist that are displaying a localized form of enhancement in the collinear limit, it 
turns out that the leading contribution is enhanced over the full comoving momentum 
triangle domain. In other words, the non-Gaussian signal is not of the enfolded type in 
this particular case and no enhancement sensitivity can be lost by the integration over all 
triangles. This dominant term to the three-point function is rapidly oscillating, causing 
the sign of the bispectrum to oscillate as well. This oscillating sign feature implies that the 
leading contribution to the scalar product with the currently available templates is severely 
suppressed and are not sensitive to the full (kirjo) 2 enhancement. Instead, subleading order 
\kir]o\ terms also contribute to the scalar product with the local or equilateral template. The 
endresult is that the local and equilateral templates only probe a linear fcir/o enhancement. 
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The details of the 2d and 3d fudge factor are complicated for the values of | Ari 770 1 = M/H of 
interest, but the numerically determined 2d fudge factor for M/H = 10 3 was used to put a 
constraint on the Bogoliubov parameter of order 10~ 2 , close to the bound derived from the 
two-point power spectrum. The main message however should be that improved templates, 
sensitive to the oscillatory nature of the dominant contribution to the bispectrum, would 
considerably tighten these constraints. The oscillatory nature of the signal in momentum 
space suggests that specific, perhaps observable, features could appear in 3d and 2d position 
space. In any case it would be worthwhile to generalize the range of available non-Gaussian 
shapes that can be compared to the data, including oscillatory signals, which we hope to 
report on in the future. Theoretically at least, for an optimal template, this would lead 
to a limit on vacuum modifications orders of magnitude stronger than the bound obtained 
from the two-point power spectrum, which would be quite remarkable. 

One important general conclusion supported by our results is that higher derivative 
corrections, which on general grounds are always expected to be present, are extremely 
sensitive to departures from the standard Bunch-Davies vacuum state. Throughout this 
paper we assumed that the combination | A;i 770 1 is independent of the actual comoving mo- 
menta involved and equal to M/H, in the spirit of the New Physics Hypersurface approach 
to vacuum state modifications. The reason for this was scale invariance of the bispectrum, 
which we relied on to allow for comparison with the available (scale-invariant) template 
shapes. Fixing t]q instead, as one would do in a Boundary Effective Field Theory approach 
to vacuum state modifications, immediately results in a scale-dependent bispectrum. It 
would be interesting to study such scale-dependent scenarios, requiring more general anal- 
ysis tools [57] , and determine to what extent (future) analysis of 3d large scale structure 
or 2d CMB data can constrain bispectrum departures from scale-invariance. As reported, 
the bispectrum or three-point function is extremely sensitive to initial-state modifications 
in the presence of a higher derivative operator, and there is no reason to think this could 
not similarly be true for all higher n-point functions. A more general perturbative analysis, 
including higher n-point functions, might lead to a hierarchy of (theoretical) constraints 
on vacuum state modifications, perhaps pointing to the standard Bunch-Davies state as 
the only consistent possibility in practice. We hope that ongoing future work in this di- 
rection can further help us understand and identify the phenomenological and theoretical 
constraints on the vacuum state ambiguity. 
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A. Enfolded triangles 



We have shown in section 2 that non-Gaussianities due to initial-state modifications for a 
canonical single inflaton action enhances in the collinear limit, corresponding to enfolded 
or squashed triangles. Though the precise shape does not overlap perfectly, here we will 
translate our proposed template to an estimator to measure the associated amplitude /§l • 
The estimators for different shapes of non-Gaussianity can be written as follows [37] 

/NL Cal = (i^nSi + (F- 1 ) 12 5 2 + (F- 1 )^ + (F- 1 ) 14 5 4 
f*f = (F- l ) 22 S 2 + (F-^tSi + (F-% 3 S 3 + (F- 1 ) 24 5 4 
/nl = (F-%3S 3 + {F~ 1 )ziSi + (F-% 2 S 2 + (F- 1 ) 34 S' 4 
bsrc = (i ? ~ 1 )44^4 + (F- 1 )^ + (F-% 2 S 2 + (F~% 3 S 3 . 

Here Fij represents the Fisher matrix and is inversely proportional to the covariance, 
the overlap, between two bispectra. In case one has a Gaussian likelihood and only one 
parameter to fit, the Fisher matrix is equal to the inverse variance; F^ = l/cr^, with a the 
fitting parameter. Here it is given by 

p, . _ \ * hhh hhh 
3 / j ft ft fi 1 

2<h<l 2 <l 3 ^h^h^h 

which is practically equal to eq. (6.8). Once more, the -B}*] 2 ; 3 are the theoretical bispectra 
of the various non-Gaussian shapes. Ci represents the total angular power spectrum, which 
contains both the CMB signal and the noise, i.e. C\ = Cf mb bf + Ni. The 6/ is the beam 
transfer function, which is instrument dependent. If the beam is Gaussian it has the form 
bi (x exp[-/ 2 ^], where a b = 0.425FWHM. 

The estimators for the point sources, the local and the equilateral are not significantly 
modified when multiple parameters are fit simultaneously [37] , and from the calculations of 
the various cosine in the paper we can assume this should hold for our template proposal, 
although to a lesser extent due to the small overlap between the local and enfolded shape. 
In the assumptions we can neglect this overlap, the Fisher matrix only has diagonal terms. 
Therefore 

^•local _ a / r? /.equil _ q / rp fcni _ q I TP f. c IT? 

JNL —^1/^11) /NL —^2/^22) JNL — "->3/^33) t>src — <-W ^44- 

Note that one can now directly compute the pre- factor oc F^ 1 , without first inverting the 
full Fisher matrix. 

In this paper we have proposed a template for the enfolded or squashed triangles, that 
should measure 



F(h,k 2 ,k 3 ) = 6A| 



1 3 
+ 2 perm + -™ 

12 1 2 3 

, A , 3 + 5 perm 

r\i\ f\i<2 /vg 



(A.l) 
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Note that if we took a slightly different template, e.g., choosing a 4 instead of a 3 in the 
equation above, the following will still hold, and one simply needs to replace this 3 with a 
4 from now on. 

To quantify to what extent the enfolded template gives complementary information 
once applied to the data, as compared to local and equilateral templates, we should consider 
the scalar product between the different templates. As pointed out in the main text, the 
starting point for deriving the most optimal enfolded template shape function is F cni = 
—F eqml + cjk\k\k\, in terms of a general parameter c. We will plot the cosine of this 
template distribution with the equilateral template as a function of c, with c running from 
to 4. The result is shown in figure 7. If we demand the template to have a definite (positive) 
sign, one should really only consider c > 1. In that case it should be clear that the optimal 
value, i.e. the smallest cosine equal to Cos(i ?enf , _F ec i ml ) = 0.49 , is achieved for c = 1, as 
claimed in section 5. For the cosine with the local template we find that it is more or less 
independent of c, as exemplified by the fact that Cos(F cnf , F local ) c=3 /Cos(F cnf , F local ) c= i ~ 
1.06 and growing ever slower for larger c. For the cosine with the local template we find, 
for c = 1, that Cos(i ?cnf , F local ) = 0.67. This is quite large, which we should have expected 
since we can imagine the local template to be a special case of the factorized enfolded 
distribution, for which only the endpoints of the line %2 + X3 = 1 are maximal, versus 
the whole line for the enfolded template. For completeness let us also mention the cosine 
between the local and equilateral template Cos(F equi1 , F local ) = 0.41. 

Any deviations from scale invariance can simply be taken into account by replacing 
power of n with n — (n/3)(n s — 1), with n s the spectral index. If one forgets about the 
pre-factors and divides out the k\ dependence, the shape can plotted as a function of k^/ki 
and k^jk\. This is shown in figure 4. 

Indeed the template shape maximizes when k\ = k2 + k%. However it does not blow up 
at this limit, which is the case when the denominator is proportional to k\ — k2 — k$. Such 
behavior can (possibly) not be mimicked, using factorizable templates, Therefore such a 
denominator would not be allowed, since a function with such a denominator can not be 
split up into functions of individual comoving momenta k\, k% and k^. 

There exists another approach [58], in which such a denominator is written as follows 



Now one has an integral over exponentials, which are factorizable. For eq. (3.2) this would 
imply 



Obviously this introduces another integration, increasing computational time one is trying 
to win by factorizing. In [58] it was shown that the double integration can be done rather 
quickly for an equilateral template. In our case, this might not be possible because of the 
large number of oscillations in the cosine, which we expect to require a large number of 
quadrature points, when one replaces the integral over t by a weighted sum. For now, we 
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Figure 7: The plot shows how the Cos(F enf , F cc J uil ) chan ges as a function of c. The smaller the 
value of the cosine, the more distinct the two shapes are. As expected for c = the cosine is -1, i.e. 
F onf = _^equii^ Note that f or < c < 1 the value of the cosine crosses zero. However these values 
of c can not be used, because the sign of the three-point function should be definite. For < c < 1 
this is not the case and the reason for getting a smaller value for the cosine is due to cancellations 
between positive and negative parts of the shape function. 

will focus on the enfolded template shape of eq. (A.l) and leave the investigation of the 
method above for future work. 

One can define the following maps 



Here Aj(fc) is the photon transfer function, introduced in section 6, which is used to com- 
pute the (theoretical) angular power spectrum C\ = (2/tt) f k 2 dkP$(k)Af. P^k) is the 
primordial power spectrum P$(fc) oc k ris ~ l /k 3 . It can be seen that all maps have a differ- 
ent primordial power spectrum dependence (and consequently different dimension). These 
maps are required to set up an estimator that has the same k dependence as the template 




(A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.l). 
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Using eq. (A. 3) through (A. 6) one can construct 4 'filtered' maps (recall that AT(ra) = 

E«m aimXlm{n)) 

A(h,r) = a ii r )j; a irnYi m {n), (A.7) 

1=2 m=-l 1 

Imax I i 

B(h,r) = V V Pi{r)+a lm Y lm {n), (A.8) 



1=2 m=-l 

Imax I 



C(h,r) = Y,Y1 li{r)j;a lm Y lm {h), (A.9) 

1=2 m=-l 1 

Imax I i 

D{n,r) = J J S x {T)-la lm Y lm {h). (A.IO) 

1=2 m=-l 1 

The sum over I runs from 2 to l ma ,x> since the monopole and the dipole are hard /impossible 
to measure and WMAP (or any other instrument for that matter) can only measure up 
to a certain l m ax based on the instrument's technical limitations. Now we need to set up 
a bispectrum that has the 'same' comoving momentum dependence as the template. For 
reasons that will become clear later, it is convenient to define the bispectrum related to 
local shape (/nl* 1 )- The local shape is local in real space and its template is exact, that 
is, the theoretical shape is equivalent to a factorized template 

F(h,k 2 ,k 3 ) = fhT 1 ^ + TO + • ( A - U ) 

It can be seen from the template that the shape is proportional to a product of power 
spectra (n s = 1). We can thus use cyclic product of the angular maps (A. 3) and (A. 4). 
The (angular averaged) bispectrum can be written as 

/•oo 

B hhl 3 = 2I hhh / r 2 dr K (r)A 2 (r) h (r) + cycl. perm] . (A.12) 
J o 

One integrates over the comoving distance. The sampling rate (in r space) depends on the 
behavior of the transfer function A;. In addition, Ii t i 2 i 3 is known as the Gaunt factor and 
is a product of the solid angle integration of (angular averaged bispectrum remember) the 
Y im . It is given by (eq. (6.4)) 



^^ r i+i)( \ +i)(2i3+i) (oooj- < a - i3 » 

The first term in the enfolded shape is equivalent to the local shape, so we can use Bj°f^ 
to express (partly) Bf^^ . The rest is obtained via carefully combining products of the 
angular maps (A. 3) through (A. 6) 

B hL = 3S£wa + &I hhh [ r 2 dr [- ((3 h (rh l2 (r)5 h (r) + cycl. perm) + 36 h (r)6 h (r)S l3 (r)] 
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Now one can easily set up the skewness estimator as explained in [37,44,46] 

cenf c i olinear 

*3 — *Jprim ~r *Jp r im ' 

where the first term is simply the term that represents the shape of the bispectrum (the 
cubic term, cubic in the filtered maps), while the second, the linear term (linear in the 
filtered maps), is added to minimize the effect caused by the inhomogeneous noise that 
breaks rotational invariance (Ci will have off-diagonal terms). The linear term should be 
constructed such that it minimizes the variance of the estimator. If the linear term is 
constructed as follows, this can indeed by achieved. One first has to derive the filtered 
maps A, B, C and D of eq. (A. 7), (A. 8), (A. 9) and (A. 10) that can be used to set up 
the cubic statistic estimator, 5 pi -i m . Subsequently one takes the Monte Carlo average, 
(S'prim)MC- Now let us suppose that S pr i m is constructed out of 3 filtered maps A, B and 
C. One can apply Wick's theorem to rewrite the average of a cubic product (ABC)mc = 
{A)mc(BC)mc + (-B)mc(^C)mc + (C)mc(AB)mc- Next, remove the MC average from 
the single maps and replace the maps within brackets with simulated maps. Our linear 
estimator becomes: ^(5 sim C s i m ) M c + £(^simC si m)MC + C(A sim B sim ) M c- If we apply this 
trick and apply weighting functions we get 

5 local = 4vr y r 2 dr J ^ (A(n, r)B 2 (h, r) 

-2B(h,r){A sim (h,r)B sim (h,r)) MC - A(h,r){B 2 im (h,r)) M c) , (A.14) 
for the local estimator and 

genf _ 3^1°°^ 

2 , f d 2 n 



+24vr / r z dr / [(-B(h, r)C(h, r)D{h, r) + B{h, r)(C sim (n, r)D sim {n, r)) MC 

J J w 3 

+C(n, r) (B sim (h, r)D sim (h, r)) MC + D(n, r) {B sim (n, r)C sim (n, r)) MC ) 

+ (D 3 (h,r) - AD(h,r)(D 2 sim (h,r)) MC )] (A.15) 

for the enfolded estimator. Here W3 is sum of the weighting functions cubed 

w 3 = [ d 2 nW 3 {h). 



The cube is a result of the fact that one has 3 a; m 's in each cubic product of the filtered 
maps in the skewness estimators. In real space a mask is simply a multiplication (i.e. one 
can multiply each AT(ra) with either zero or one). However, this becomes a convolution 
in Fourier space. Consequently we have an integral over the solid angle d 2 n. If there is 
uniform weighting (that is, each pixel is masked or not) W{n) = M(n), the mask function 
and the integral becomes 

■^3 = ^fsky 

with f s ky the sky cut (in fact it is the opposite, it represents the percentage of sky that 
remains after masking). In the latest WMAP analysis they have not used a uniform 
weighting, but a "combination signal-plus-noise weight", which turns out to be optimal 
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for the analysis of equilateral shaped bispectra, while the local shape is barely affected by 
simply using the uniform weight. It should be checked to what extent uniform weighting 
changes the estimates of the enfolded shape, compared to the more advanced combined 
weighting used by WMAP team [37] . 



B. Modified initial-state bispectra 

We start with eq. (3.11) in [54], from which we can extract the interaction Hamiltonian 
from the canonical effective action minimally coupled to gravity 

/ . \ 4 



V 



Hl = ~W d xa{v) H C 9 C ■ 



As carefully explained in [54] the three-point correlation function (CfciCfoCfcs) can be writ- 
ten as (to first order in the interaction Hamiltonian Hj) an integral over the free field 
correlator { Cfci Cfc 2 Cfc 3 -f^/ C 1 ? ) ) ) where the free field correlator can be expanded in a product 
of two point correlators (i.e., Wightman functions) via Wick's theorem. Consequently it is 
straightforward to show that the three-point correlation function in the case we consider 
the interaction Hamiltonian of eq. (B.l) is given by 

(CkMks) = - i ( 27r ) V3) (E^)]^(l : ) x (R2) 

o i 

dr/a^r,)-^^ (O,77)0„G£ 2 (O, V)9 V G> 3 (0,7/) + perm + c.c, (B.3) 

as can be found in [54]. Here the Wightman functions are defined as followed 

(CfcCte) = (2^)V 3 )(fc! + k 2 )G>(r,, V '). (B.4) 

The Wightman functions can be found by solving the e.o.m. of the inflaton field minimally 
coupled to gravity 

GlM) = f^^ + ^X 1 -ik V ')e- ik ^'\ (B.5) 



Consequently we compute 



G>(0,f,) = ~^(l-ikr,)e^, (B.6) 

and 

d ri G>(0,r ] ) = - I !^^-^-e^, (B.7) 
^2 2ka(r]) 

with a{rj) = 1/r/H during inflation in the assumption H ~ 0. 

Next we can express {CkiCk 2 Ck 3 ) m terms of these Wightman functions. However, we 
are interested in the case where we are not in the BD vacuum state. To first order in the 
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Bogoliubov parameter (3k, what happens is that one of the solutions to the equation of 
motion 'picks up' a minus sign in k. It is easy to incorporate this minus sign to find 



= -it*?** (Eft) 4rwH/^^ e " J '' +c - c ' 

(B.8) 

Note the factor of 6 is the result of the 6 possible permutations, while the sum is the 
result from implementing an initial-state modification for all of these 6 permutations, kj = 
kt — 2k j, with kt = k± + k2 + k%. The integral can be easily performed as well as adding 
the complex conjugate part. In steps 

• f° * ik v 2(cos(Lry ) - 1) 
-i x / drje™^ + c.c. = — \JJ^L L 



rio kj 



resulting in 



(B.9) 

This is the result we have used in eq. (3.1) of section 2. In case we assume that the 
enhancement occurs when kj — > 0, we can apply this limit to the above expression to find 

= -wv> (Eft) w P ukf)$^ 3A p ne{ ' %,) 



x-^„§ + 0(fc?), 



which is slightly different from the result found in [54] since they did not consider the limit 
of x — > in cos(x)/x correctly. Note that when kj = this whole expression actually 
vanishes, but does have a maximum nearby (i.e., k m ax ~ ?7o )• 
Next we consider higher-order terms in the action of the form 

£/ = v^^av^)^)) 2 . (b.io) 

It is not hard to compute the corresponding interaction Hamiltonian up to third order in 
the curvature field Q (and ( ~ — (H/c/>)54>, where = 4>{if) + 5<f)(r],x)). It can be shown 
that 

Hi = ~ J d*xa( v ) £ (C 2 - (9,C) 2 ) • (B.H) 

In similar fashion we can write down the three-point correlator in terms of integrals over 
products of Wightman functions, while Fourier transforming to k space 

(C^CO™ = -*(2vr)V 3 ) ^ (|)7 d V a(r,) 

(8^1(0,^01(0^)8,01(0^) + h ■ fc 2 G> (0,^)01(0^)8,01 (0,7?) + perm) , 
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and its complex conjugate. Here the dot product comes from the Fourier transform of two 
partial spatial derivatives in the interaction Hamiltonian of eq. (B.ll). This is a much 
'larger' three-point correlator, so let us compute it in such a way that we do not lose track 
of all different components. The best approach is to first compute all the different terms in 
the integral and then group these in proportionality to rj (i.e. oc rj , r\ and rj 1 ) . In addition 
we can compute the whole 'pre-factor' independently. It is not hard to show that it is given 



(B.12) 



and therefore 





no 



+kj ■ kj + ikj +2 (l + ikjrj)(l - ikj+irj) 
+kj ■ k j+2 kj +1 (l + ikjT})(l - ik j+2 r]) 
+k j+ i ■ k j+2 k](l - ik j+1 r])(l - ik j+2 r]) 



(B.13) 



where the j's are cyclic in 1, 2 and 3. 

Let us first rewrite the dot product using the triangle vector constraint k\ + k 2 + k% 
Using this constraint we can deduce that 

1 



0. 



2^1+2 



ki 



kj+i)- 



If we write the integral as J*^ dr/ Y2j P% S{k\, k 2 , k$, r/)e lfcj ' r? we obtain 



1 



S(ki,k 2 , k 3 , rj) = +- [(fe| +2 - k] - k 



+ijk]+2 + kj + i 



kj+2)kj+i 



+{kj — kj + i — kj +2 )kj] + — [(^|-|-2 ~~ kj — kj +1 )kj +2 (kj — fcj+i) 

~ kj — kj +2 )kj +1 (kj — kj +2 ) 
-ik) - k] +1 - kj +2 )k](k j+1 + k j+2 )] 



+y [kjk j+ ik j+2 ((k j+2 (k] +2 - k) - fc? +1 ) + k j+ i(k] +1 - k) - k 2 j+2 ) 



—kj(kj — kj + i — kj +2 )) + 6kfk%kl] . 



(B.14) 



This can be rewritten as 

1 

— i 
2 



S(ki,k 2 ,k 3 ,rj) 



h Y\ ki + — kj (2kj + ik j+2 {kj+i + kj +2 )' 



-k t (kj(k] +1 + kj +2 ) + (kj+i + k j+2 )(2k] +1 + 3k j+1 k j+2 + 2k] +2 ) 
-kj(k j+ i + k j+2 )k^j 

-Akj + ik j+2 {kj + i + k j+2 ){k 2 j+1 + k] +2 + k j+1 k j+2 )] 



' 2 



h Y\ ki{k 2 - Akj + ik j+2 



(B.15) 
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which is very similar to eq. (3.30) in [54] except for the term linear in 77, for which we 
find additional terms. Since the integral runs over conformal time rj, all that is left is 
to compute the different integrals in terms of rj. Again, it is useful to pre-compute the 
following integrals 



x 



/" 

J vo 



drj e lk jV _|_ c c 



-i x 



dr}{ir\)e %k ^ + c.c. 



2(cos(fcj?7o) - 1) 



2r/ sin(A; i 7?o) . 2(1 - cos(%j/o)) 



+ 



</0 



-l X 



d-qrfe^i 11 + c.c. 



2t?§cos(A;j7/o) 4r] sm(kjT] ) 4(1 - cos(kjT]o)) 



no 



+ 



Since we now have terms that are proportional to kj one would expect some terms to 
diverge in case kj — * 0. However, this proportionality only appears in the last term of 
eq. (B.15), which is proportional to kj on itself. Consequently we lose a factor of kj after 
multiplication, just enough to make that term finite for kj — > 0. Similar argumentation can 
be applied to the other terms, once we realize that the limits of sin(x) jx and (1 — cos(x))/x 2 
are finite in the limit x — > 0. The three-point correlator therefore becomes 



(CfeiCfe2Cfc 3 )nBD 



(1 - cos(fcj?7o)) 



T] s'm(kjr] ) 



-h Y[k + Ai(k) + i (A 2 (fc) + 2B(k)) 



kjAi + A 2 (k) + 2B(k)) + v $cos(kjr) )B(k) 



(B.16) 



where 
Ai(kj, k 



'j+i, kj+2 



) = 2k j+1 kj +2 (kj+i + fej+2)' 



-4-2(fej)%+l)%+2) 

B(kj,kj + i,kj +2 ) 



+(/c j+ i + fcj+2)(2fc| +1 + 3/cj + i/cj +2 + 2fcJ +2 )] - kj{kj +1 + k j+2 )k 2 
-4A; j+ i/c j+ 2(fcj + i + kj +2 ){k 2 j +1 + fc| +2 + k j+1 k j+2 ) 
\h(k 2 t -Akj + ik j+2 ). 



Unlike in the simple case, and aside the comoving dependence of the nominator, this three- 
point correlator has terms that are proportional to (1 — cos(kjr]o))/kj and s'm(kjrjo)/kj. 
Therefore we can take the exact limit kj = 0, and have a no n- vanishing result. The easiest 
way to see what happens is to look at the integrals. In the limit kj = the first integral 



- 36 - 



and the last integral vanish. Consequently, we can write 



rf 

x-£-A2(kj,k j+1 ,k j+2 ) 

„2 

x y (-4fc j+ i/c i+ 2(A; i+ i + %+ 2 )(fcf + i + fc| +2 + k j+1 k j+2 )) , 

which is equivalent, up to a minus sign, to what the authors of [54] found. However, we have 
shown here and have argued in section 5 that this limit does not represent a maximum. 
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